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Abstract 

We obtain rigorous upper and lower bounds for the error in the recent approxima- 
tion for 7r proposed by Chakrabarti & Hudson. 

1 Introduction 



In 2003, Chakrabarti and Hudson pQ proposed a new formula in the spirit of Archi- 
medes for the computation of tt, namely 

OO- i 

tt « — {327T n + 4II 2n - 6a n }, (1.1) 

where ir n is the perimeter of an inscribed regular polygon of n sides in a circle of perimeter 
tt, U n is the perimeter of a circumscribed regular polygon of n sides in a circle of perimeter 
tt, and a n is the area of an inscribed regular polygon of n sides in a circle of area tt. 

Archimedes, in The Measurement of the Circle (MC), see [2], obtains his own approxi- 
^ 1 mation by starting with the (evident) inequality 

7T n < 7T < Il n , (1.2) 

and takes n = 6,12,24,48,96, respectively. He performs a brilliant tour de force of ma- 
nipulating and rounding the inequalities resulting from continued fraction expansions of the 
square roots which arise in the computations of 7T n and n n , whence he obtains his justly 
famous bounds 

10 1 , . 

3— <tt<3-. 1.3 
71 7 V 1 

It is generally recognized, see [6], that the extant MC is a post- Archimedean revision 
of Archimedes' original and far more comprehensive treatment of the circle. For example, 
Heron, Metrica I, 32 [3], quotes a theorem of Archimedes from (MC), which is not present 
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in the extant version. The theorem states that the area of a circular sector exceeds four-thirds 
the area of the greatest inscribable triangle. In terms of the notation above, this inequality 
affirms that: 

7T > i(47T 2 „-7r n ). (1.4) 

The right hand side of fll.4j) is a convex combination of 7r 2n and 7r n since 3 — | = 1. More- 
over, its form argues for the fact that the original treatise of Archimedes studied convergence- 
improvement inequalities derived by means of geometric theory. Unfortunately, such the- 
oretical aspects apparently proved too subtle for the uses of later commentators such as 
Heron. 

Observe that the formula of Chakrabarti & Hudson (11. ip . too, is a convex combina- 
tion of its terms. 

Such convex combinations did not reappear until the sixteenth and seventeenth centuries 
in the work of Snell [7] and Huygens [5j. The former stated and the latter proved the 
following upper bound inequality: 

^<§7r„ + in n , (1.5) 

which is a beautiful generalization of Archimedes' original inequality. Both Snell and Huy- 
gens use it in the form 

7T W §7T n + |n n (1.6) 

and they justifiably extol the convergence-rate improvement produced by this convex com- 
binations of quantities already computed. 

Now, Chakrabarti and Hudson ask what happens if we use the convex combination 

7r^|a n + |A„ (1.7) 

to approximate n, where A n is the area of a circumscribed regular polygon of n sides in a 
circle of area tt. They prove the very interesting result that 

lim f n + f; - 71 = I , (1.8) 
~a n + i A n - it 8 

which shows that the areas converge much more slowly than the perimeters to it. 

Then they suppress the limit notation, treat (11.81) as an equation for n and when they 
reduce it algebraically they obtain the approximative formula (11. ip . Passing to trigonometric 
functions, ( II. ip becomes: 

^s{ 32sin © +4tan GD - 3sin ©} sn(,i) ' (L9) 

Finally, the authors offer some numerical studies of the accuracy of ( 11.91) . 
The investigation of the authors needs to be completed in several areas. 

• They do not state whether the approximation (II .9p is in excess or in defect. . . in fact, 
we will prove that it is in excess. 
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• They do not develop any rigorous error analysis. . . no upper bounds for the error nor 
lower bounds; we will present such bounds. 

• They offer a non-standard definition of accuracy, which they call "precision". This 
makes it difficult to compare their numerical results with standard error studies. We 
will present the standard definition of "correct significant digits" and reformulate the 
discussion of the accuracy of (11.91) in light of our error bounds. 

2 Error Analysis 

We will prove the following theorem: 

Theorem 1. If n ^ 32 then the following inequality is valid: 



1 f n \ 6 < n(n) — 7T 1 fit 



105 Vn/ 7r " 104^ \n 



(2.1) 



where the lower-bound constant is the best possible. 

105 

An immediate consequence is: 

Corollary 1. The approximation tc pa U(n) is in excess. □ 

The inequality (12. ip bounds the relative error in the approximation IT(n) ~ re. It is well 
known that an approximation has n correct significant digits iff the relative error does not 
(!) 

exceed jg^' see El- Therefore, we can say: 

Corollary 2. The approximation n m U(n) has about (61og 10 n — 1.27) correct significant 
digits. □ 

Now we turn to the proof of (12. II) . We will use the MacLaurin expansions of the 
functions involved. Define: 

f(x) : = — (32 sin a; + 4 tan x -3sin2x). (2.2) 

Then we see that 

n(») = /(£). (2.3) 

Since the sum of a convergent alternating series is bracketed by two consecutive partial 
sums, we obtain: 

Lemma 1. The following inequality is valid for all x ^ 0: 

357\ /3579 
X> Xj Xj \ I Xj Xj Xj Xj 



^ '-3! + 5! -71) <32sinx<32^-- + - -- + -). (2.4) 

□ 
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The MacLaurin expansion of the tangent function is not an alternating series. But we 
can still use the Lagrange form of the remainder to obtain: 

Lemma 2. The following inequality is valid for ^ x ^ ^: 



32 ' 

J 1 3 2 5 17 7 62 ,\ A 

\ X+ 3 X + 15 X + 3lT + 2835 X ) <4tanX 

J 1 3 2 5 17 7 62 9 1 lA 

<4 x + -x 3 H x 5 H x 7 H x 9 H x 11 ). (2.5) 

V 3 15 315 2835 85 / v ' 

Proof. The MacLaurin expansion of order 11 of the tangent function is 

X 2 17 62 X 

tanx = x + -x 3 + —x 5 + — x 7 + —x 9 + —R(9 n )x n , (2.6) 



where 



d 



n 



R(x) := — ^(tanx) = 256(tan 2 x + 1)(155925 tan 10 x + 467775 tan 8 x + 509355 tan 6 x 
ax 11 

,4 _ , mine j.„_2 



+ 238425 tan 4 x + 42306 tan^ x + 1382) 
and ^ On ^ ^. The function i?(x) monotonically increases in the interval, whence 

353792 = i?(0) < R(6 n ) < R (—) = 469223.9941 . . . (2.7) 

V 32 / 

and we conclude (since x ^ 0) that 

^ — :R(6 n )x n ^ — ■ 469223.9941 . . .x 11 = x 11 < — x 11 . (2.8) 

11! v ; 11! 85.06. . .85 y J 

This completes the proof. □ 

Finally, as in the inequality (12.4ft we obtain: 
Lemma 3. The following inequality is valid for all x ^ 0; 

f (2x) 3 (2x) 5 (2x) 7 (2x) 9 l 
- 3^ 2x - + - + } < -3 sin 2x 

\ 3! 5! 7! 9! J 

f (2x) 3 (2x) 5 (2x) 7 (2x) 9 (2x) n l , . 

□ 

Now that we have set up the technical inequalities necessary in our main proof, we enter 
into its details. 
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Proof of the main theorem. If we substitute the inequalities (12. 4p . (12. 5 p and (12. 9p into the 
formula for f(x), (12. 2p . we obtain 

1 6 1 8 2776 10 f( X ) -, 1 6 1 8 10531 10 

-x 6 + x 8 + —x 10 < - 1 < x 6 + x 8 + x 10 . (2.10) 



105 360 2338875 n 105 360 26507250 

Now, the left-hand side 

1 6 1 8 2776 10 1 ef- 7 2 2776 A 1 6 /n^ 

x 6 + x 8 + x 10 = x 6 <^ 1 + —x 2 + x 4 } > x 6 (2.11 

105 360 2338875 105 \ 24 22275 J 105 v ; 

for all positive x. So, we have proven the lower bound (12 .ip and that the constant 



cannot be replaced by a bigger one, i.e., it is the best possible. 
For the upper bound, we observe that for ^ x ^ ^, 

x 6 + —x 8 + 10531 = —x 6 \ 1 + —x 2 + 10531 x 4 



105 



105 360 26507250 101 | 2 1 252 150 

< 



105 \ 24 V 32/ 252450 V 32/ 



1 6 



104.705... 104 1Q 

which is the upper bound presented in (12.11) . 

This, together with f ]2. 31) complete the proof of the theorem. □ 

3 Numerical Error Studies 

Chakrabarti and Hudson present a table of numerical studies of the error in their approxi- 
mative formula (II. 8p . They define: a number a has precision n if 

1 / s 

\a-7T \< . 3.1 

1 1 10™ K J 

We already pointed out the standard definition (see [1]): an approximation N approximates 
the true value N with n correct significant digits if the positive relative error 



liV-M 



N <w- (3 - 2) 

Of course the two definitions will coincide sometimes, and sometimes not. However, we note 
that 

1 la-d 1 (i) 

Qt — 7T < 1 < < 1AL 

1 1 10" 7T 7TlO n 10™ 

so that an approximation with precision n is always correct to n significant digits. It is 
sufficient, but not necessary, since if 

1 let — 7r| (tt) 
< < ' 



7TlO n 7T 10 n 

then a will have precision n — 1 but still have n correct significant digits. 
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4 Earlier approximations 



We can rewrite the approximation 1 



/(*) 



7T 



sin x 



x 



as follows: 



15 cos x 



2 + 16 cos x — 3 cos 2 x ' 



(4.1) 



x 

where we know that the approximate value is smaller than the true value by about ^g^ - 
sin a; 



Now, 



x 



has the following continued fraction expansion: 



sin x 

x 



i|sin 2 f 



M sin2 f 



(4.2) 



1 - 



|| sin 2 f 



1 - 



34 • 2 x 
7-9 blil 2 
1 



The first few convergents are: 



Pi 
Qi 

El 

n 



2 + cos x 



9 + 6 cos x 



14 + cos x ' 

51 + 48 cos x + 6 cos 2 x 
80 + 25 cos x 



error 



error 



error 



180 
1 

2100 
1 

"44100 



x 4 + - 



x° + 



X* + 



(Snell) 
(Newton) 



(4.3) 



where we have used sin 2 | = ^(1 — cosx). (For all these results, see [Sj.) Every one of these 
approximations is larger than the true value and is the best possible. 

The Chakrabarti-Hudson approximation (14. ip has an error about twenty times greater 
than that of Newton's formula, although both are of order o(x 6 ), nor is not as simple. Indeed, 
it has the formal complexity of the third convergent without the latter's extraordinary accu- 
racy. Nevertheless, it is the only approximation in defect, its accuracy is still quite good, and 
it is interesting that such an ad hoc derivation produced such an intriguing approximation. 
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